New Universality Class at the Superconductor— Insulator Transition 
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We study dynamic properties of thin films near the superconductor - insulator transition. We 
formulate the problem in a phase representation. The key new feature of our model is the assumption 
of a local ohmic dissipative mechanism. Coarse graining leads to a Ginzburg-Landau description, 
with non-ohmic dynamics for the order parameter. For strong enough damping a new universality 
class is observed. It is characterized by a non-universal d.c. conductivity, and a damping dependent 
dynamical critical exponent. The formulation also provides a description of the magnetic field-tuned 
transition. Several microscopic mechanisms are proposed as the origin of the dissipation. 



PACS numbers: 74.76.-w, 74.50. -f-r 

Granular superconductors and Josephson junction ar- 
rays behave similarly in many respects. This is because 
the key degrees of freedom are thought to be exclu- 
sively bosonic in nature, the Cooper pairs of the un- 
derlying electronic problem. Quantum phase transitions 
are present in the ordered arrays, when superconductiv- 
ity gives way to a gapped (Mott-) insulator. This is a 
direct consequence of the uncertainity relation between 
phase and charge degrees of freedom . When disorder 
is present, an additional insulating phase may appear, ex- 
hibiting glassy behaviour In two dimensions, early 
experimental studies reported the conductivity a, to as- 
sume a universal value A theoretical explanation of 
this universality was proposed based on scaling consider- 
ations j^]. However, subsequent measurements observed 
critical resistivities, differing as much as tenfold 
The main motivation of our work is to develop an under- 
standing of this apparent non-universality. 

In this paper we develop a Ginzburg-Landau- Wilson 
(GLW) formulation for the superconducting order- 
parameter. We emphasize the importance of electronic 
excitations to account for the experiments. Our main re- 
sult is the observation of a new universality class at the 
superconductor - insulator transition. The conductivity 
at criticality is found to be non-universal, and the dy- 
namical critical exponent depends on model parameters 
for strong enough dissipation. The non-ohmic dynamics 
of the order parameter leads to a power-law behaviour 
of the conductivity at low frequencies in the insulating 
phase. The magnetic field-tuned transition will be also 
discussed. 

Our starting point is an array of superconducting is- 
lands or grains connected by Josephson junctions. The 
jth gj-g^jj^ jg characterized by the phase (fi of its supercon- 
ducting order-parameter. It is also important to include 
the low energy electronic degrees of freedom. Integrating 
out the electrons leads to a dissipative dynamics for the 
local phases. (For a more detailed physical motivation 
see below.) One arrives at the Euclidean action 
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The smallness of the grains manifests itself in a charging 
energy U which is the energy cost of transferring an ex- 
tra Cooper pair on the island. The Josephson coupling 
between neighboring islands is denoted by J. If the low 
lying electronic excitations have a finite density of states, 
then the damping will be ohmic; its Fourier transform is 
thus given by a{iv) = ao\u)fj,\/2Tr. The locality of the 
damping is reflected in the fact that the dissipative term 
in Eq. (|l]) couples the phase in a single island at different 
times. Additional non-local charging or dissipative terms 
do not change qualitatively the results and for simplicity 
will be dropped. 

A Hubbard-Stratonovich transformation is employed 
to decouple the Josephson term [|l|. This introduces the 
complex order-parameter field ip such that its expectation 
value is proportional to that of exp(iip). This yields the 
GLW action which in two dimensions takes the form 
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The correlation function g is given as: g{T) ~ 
(exp{i(^i(r) — iipi{0)})o, where (• ■ •)o is an expectation 
value taken with respect to the single site Gaussian part 
of the action of Eq. (Q). In the presence of local damp- 
ing g{T) decays algebraically in time (~ ,--2/00^^ ^j^g 
Fourier transform for small frequencies reads 
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16J/[C/(2 — q;o)]. To summarize the low frequency form 
of the GLW action reads: 
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where e = 1— go- Clearly the model exhibits a non-ohmic 
dissipative dynamics, reducing to ohmic damping only in 
the special case s = 1. Surprisingly, the ohmic damping 
in the quantum phase model of Eq. (|^) yields a non- 
ohmic dynamics for the coarse-grained order-parameter. 
The quadratic frequency dependence dominates over the 
dissipative dynamics for s > 2, or equivalently ao < 2/3. 

The model exhibits a phase transition at e(ao, J, U) = 
0. The ordered phase supports long-range superconduct- 
ing order, whereas the disordered phase is an insulator. 
The phase diagram for T = is displayed in the inset 
of Fig. 1. Increasing damping shifts the phase bound- 
ary to lower values of J. The insulating phase disap- 
pears for ao > 2. In this region the dissipative processes 
completely suppress the quantum fluctuations, in anal- 
ogy with the dissipative phase transition, discussed in 
Ref. i. 




FIG. 1. Conductivity at the transition as a function of ao- 
The inset shows the corresponding phase diagram at T = 0. 
Along the solid line a is universal, whereas it is a function of 
Qfo along the dotted line. 

At r = a quantum phase transition takes place which 
is characterized by a dynamical critical exponent z. In 
the limit of weak damping, the critical behaviour is that 
of the non-dissipative models, z = 1. In the general case 
it is given by z = max(l, 2/s), and the critical behaviour 
is that of a (2 + z)-dimensional XY-model. Therefore our 
results establish the existence of a new universality class 
of the superconductor - insulator transition for s < 2. 

This quantum phase transition attracted intense inter- 
est because of the claim of a universal conductivity at 
criticality Q . To address the issue of universality we cal- 
culate the conductivity in terms of two and four point 
Green's functions [ pO[ . In the Gaussian approximation 
the four point function factorizes. An analytic continu- 



ation iuj^ Lu + iS yields the dependence on the real- 
frequency 
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where ctq = h/Ae^. The advanced and retarded Greens 
functions are given by 
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cos Y ± isign(w) sin y , (6) 



with e > 0. Without damping the frequency dependence 
exhibits a gap of the size luq = ■\/4e/C, thus this disor- 
dered phase is a Mott-insulator. With increasing damp- 
ing the Mott gap is smeared out. For s < 2 and low 
frequencies uj <^ loq ~ -^/leTC one finds 
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The conductivity shows a power-law behaviour at low 
frequency, where the power depends on the strength of 
the dissipation for s < 2. These two types of behaviours 
are displayed in Fig. 2. 
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FIG. 2. Real and imaginary parts of the conductivity as a 
function of the frequency. Without dissipation (ao = 0) the 
real part exhibits a gap (dotted lines). This gap is smeared 
in the presence of dissipation (solid lines). The damping pa- 
rameters are s = 1/2, and riu)o/e — 1. 

Of particular interest is the d.c. conductivity at the 
transition, i.e. for uq (with ujq/uj 0). As shown 
in Fig. 1 the value of the d.c. conductivity at criticality 
a* is a function of the strength of the ohmic damping for 
ao > 2/3. To summarize, the non-dissipative supercon- 
ductor - insulator transition possesses a finite basin of 
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attraction. In the region of < ao < 2/3 the dissipation 
is an irrelevant operator. It is characterized by z = 1 and 
a universal critical conductivity. For strong enough cou- 
pling to localized fermions ao > 2/3 a new universality 
class is present, with damping dependent a* and z. This 
is the central result of our paper. 

Renewed interest in the superconductor - insulator 
transition was generated by the early experimental re- 
port of a universal conductivity at criticality, with the 
value (cr*)"^ = 6.5kn Q. However, subsequent measure- 
ments on superconducting films and on Josephson 
junction arrays Q reported resistivity values distributed 
in the broad range of 2-20 kft. Our emphasis on the elec- 
tronic degrees of freedom offers a possible explanation of 
this lack of universality. 

In many experiments the transition is tuned by a mag- 
netic field B [|l^,D, which we therefore discuss next. The 
magnetic field is incorporated by introducing a vector 
potential A{x, r) in Eq. (Q). For weak frustration / ^ 1 
the lattice structure can be neglected. The magnetic field 
leads to the formation of Landau levels. The longitudinal 
conductivity is given by 
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^ , with the 
The mag- 



where Gi^^^n = [e+U!cin+l/2)+ri\uj^\'' + Ciul] 
dimensionless cyclotron frequency lJc — tt/ 
netic frustration / = $/$o is proportional to the flux 
through a unit cell. In the mean-field approximation the 
phase boundary is given by e -I- uJc/^ = 0. The fact that 
the location of the phase boundary is magnetic field de- 
pendent demonstrates, that our formulation is capable of 
capturing the field-tuned transition. 

The analytic continuation to determine the conductiv- 
ity as a function of real frequencies follows the same lines 
as in the zero field case. We explicitly evaluate the real 
part of the conductivity as shown in Fig. 3. The fre- 
quency dependence reflects the underlying Landau-level 
structure. It is smeared due to the influence of the damp- 
ing. The d.c. conducitvity at the transition is determined 
by the lowest Landau level n = 0. In the absence of dissi- 
pation it diverges. Again we find that for strong enough 
damping, ao > Ij the critical conductivity a* becomes 
finite, with the damping dependent value, as shown in 
Fig. 3. For ao ^ 1 it assumes the asymptotic form 
a* ^ l/(ao - 1). 

The previous analysis showed that the effect of the 
local damping for the phases has a great impact on the 
properties of the system. In the second part of this paper 
we discuss various scenarios which can give rise to local 
dissipation. 

The key feature giving rise to dissipation is the pres- 
ence of gapless degrees of freedom. These electronic de- 




FIG. 3. Longitudinal conductivity at the transition in a 
finite magnetic field as a function of ao. The inset shows the 
real part of the frequency dependent longitudinal conductivity 
in a magnetic field (oscillating curve), compared to the zero 
field case (smooth curve). With the parameters uic/\/^ ~ 2, 
s — 1, and r]UJo/e — 1/4. 



grees of freedom can either be of intrinsic origin, or they 
can result from a coupling to the substrate. We argue 
that both scenarios are viable candidates for the descrip- 
tion of the experimental data. 

Before reviewing the microscopic motivation we point 
out the important difference between the local dissipa- 
tion, investigated here, and the non-local forms, studied 
in Refs. ||^. When local damping is present, the number 
of Cooper pairs is not conserved: they can decay into a 
pool of normal electrons. On the contrary, if the damp- 
ing is non-local, e.g. it is induced by inter-grain electron 
transfer, then it conserves the charge within the film. Be- 
cause of this essential difference it can be expected that 
genuinely new physics arises in the present model. 

In the remainder of the paper we discuss possible ori- 
gins of the local damping. The two necessary ingredients 
for the proposed mechanism are local pools of gapless 
electrons with a temperature independent density, and 
a coupling term between these electrons and the phase 
of the order parameter, which leads to ohmic dissipa- 
tion. Once gapless electrons are present, the Andreev 
process satisfies these criteria. It describes the conversion 
of normal electrons into Cooper pairs at superconductor 
- normal (SN) interfaces. In the absence of a tunnel bar- 
rier the charge transfer across the interface is continuous. 
This results in a fully ohmic dynamics for the supercon- 
ducting phase. Therefore the proper description of this 
process is exactly the dissipative term in Eq. (|l|). For 
the case of two superconducting islands, connected by a 
normal metal, a microscopic derivation of the ohmic dy- 
namics was provided in Ref. ||l^. Adapting their deriva- 
tion for the present case of an SN structure confirms our 
statement explicitly. In the ideal case the value of the 
Andreev resistance is half that of the normal state resis- 
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tance, consequently the arising local dissipation can be 
large. We emphasize that this process involves pools of 
localized electrons which are therefore not part of the con- 
ducting path. Transport across the sample is dominated 
by tunneling between grains, resulting in a much higher 
subgap resistance. 

The remaining question is the origin of the gapless elec- 
trons. In order to give rise to ohmic dynamics, the level 
spacing of their spectrum has to be small compared to 
all other energy scales in the problem. This necessitates 
a relatively large spatial extent of their wavefunction, 
which can be estimated to be in the range of hundreds 
of Angstroms. An intrinsic scenario can be suggested 
by recalling that in the absence of superconductivity the 
system is an insulator, otherwise the experiments would 
display a superconductor - metal transition. On the other 
hand, in order to be able to support superconductivity, 
it has to be close to the metal - insulator transition. In 
this region the combined effects of disorder and enhanced 
Coulomb repulsion act as very effective mechanisms to 
break up the Cooper pairs into normal electrons In 
the vicinity of the transition the localization length of 
the electrons is strongly enhanced. Thus their spatial ex- 
tent is large, leading to a considerable density of states at 
the Fermi energy. The existence of such low energy ex- 
citations with a density comparable to the normal state 
density was observed in tunneling experiments Q . 

An extrinsic origin of the local damping mechanism 
can be due to a conducting substrate. In this case the 
unpaired electrons of the ground plane constitute the 
reservoir. It is well established that in proximity coupled 
junction arrays the metallic substrate indeed is a source 
of localized dissipation |16-18|. Superconducting films 
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are commonly grown on disordered semiconductor sub- 
strates. As was suggested in Ref. ^ the semiconductor 
substrate might be able to support metallic transport on 
short length scales. It is possible that the deposited metal 
dopes the amorphous semiconductor substrate, giving 
rise to a temperature independent density of states at 
the Fermi level. An additional remarkable possibility was 
raised in connection to the study of Ag/Ge, where the 
metal - semiconductor interface itself seems to support 
metallic transport |l9[ . 

In sum we developed a modified understanding of the 
superconductor - insulator transition. We emphasized 
the importance of localized gapless electronic excitations. 
They give rise to a local ohmic dissipation for the phase of 
the superconducting order parameter. A new universal- 
ity class of the transition was identified for strong enough 
dissipation. It is characterized by a non-universal value 
of the critical conductivity and a damping dependent dy- 
namical critical exponent. The dependence of the above 
results on a magnetic field was also analyzed. We re- 
viewed several possible microscopic origins of the gapless 
electrons, and identified the Andreev process as a likely 
source of the ohmic dissipation. 
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